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UNCLASSIHED 



ABSTRACT 



Determination of an underwater target's position using passive acoustic sensors is of 
considerable use for the Navy, both for anti-submarine warfare (ASW) and underwater 
surveillance. This thesis proposes and develops localization algorithms capable of 
passively determining the location of a transient source given some broad constraints. In 
particular, this thesis investigates the effect of the source signal uncertainty on localizer 
performance. The localization process consists of two parts. First, a time domain 
propagation modeling code determines the impulse response of the environment from all 
possible source locations to a single hydrophone. This program predicts the signal as it 
would appear at the receiver from a grid of possible source locations. Second, source 
localization results from finding the maximum correlation between the positionally 
dependent, numerically modeled signals and the actual received signal. The position of the 
maximum cross correlation reveals an estimate of source position. Using model to model 
correlation, this technique successfully localized acoustic sources in both Monterey Bay 
and Barents Sea scenarios. 
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INTRODUCTION 



The localization of transient acoustic signals is of potential use to the United 
States Navy both for antisubmarine warfare and underwater surveillance. The ability to 
localize a signal with only the received signal data and with little or no a priori 
information about the original source would be a significant advance in both these areas. 

The classification of these signals is of similarly high interest. It would appear that 
all signals of a short duration nature could have a myriad of frequency characteristics. To 
conduct a true matched signal cross-correlation, knowledge of these signal parameters 
would help immensely. Computer simulated localization experiments have shown that 
near perfect correlation was possible using the signal transmitted through the model from 
the source then correlated with a like signal sent back through the model from the receiver 
to the source location. As part of the development of this technique, this "matched-signal" 
algorithm was tested for several varied cases to confirm reliability. 

In the history of broadband source localization, most previous work involves 
methods that operate in the time domain. The earliest matched-signal work, which 
ironically came as a precursor to the recent research in matched-field theory, was 
conducted by Parvulescu* and Parvulescu and Clay,^ in which an impulse was transmitted 
by a source in a given environment, the resulting time series was recorded at a receiver 
and the time-reversed time series was re-transmitted through the source. Both laboratory 
and ocean experiments indicated that the received signal reached a peak when the source 
and the receiver were at the original positions, rather than any other positions. The state of 
computer technology in the period of these experiments limited the tests to primarily 
empirical analog efforts with no practical tactical applications. 
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More recently, Clay has re-examined the matched-signal concept as a time-domain 
matched-field method. In this algorithm, a computer propagation model predicts the 
impulse response for a number of source positions, counter to physically measuring the 
impulse response of the laboratory or ocean environment. The time-domain signal 
matching method then localizes the source. Clay also extended the method to multiple 

receivers by cross-correlating the matched-field outputs for pairs of receivers The \ 

\ 

research continued with Li and Clay applying the methods to experimental data measured 
in rigid- walled wedges and flat waveguides. Hodgkiss and Brienzo successfully applied 
the method to short-range localization of impulsive signals in deep water. Later, Frazer 
and Pecholcs proposed generalizations to the matched-signal algorithm for a single 
hydrophone by considering different norms.* Most recently, Westwood developed a 
broadband matched field algorithm that operates in the frequency domain without a priori 
source knowledge.^ 

A localizer based on the correlation of a signal with that of a signal sent through a 
propagation model, i. e., matched-signal, requires a prototype for a signal. Given that the 
transmitted transient is assumed to be an impulse in time, a impulse response function 
model of the signal through the medium would possibly suffice. As did Westwood,^ this 
thesis will show that, indeed, under given constraints of both the environment and the 
signal to be localized, that this algorithm will work. 
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II. MODELING THE OCEAN TRANSFER FUNCTION 



A. PARABOLIC EQUATION MODEL 

Under the assumptions of linearity and time invariance,* an acoustic signal at a 
receiver can be found from the convolution of the source pressure signal with the ocean 
impulse response, 

OC 

p(r^,z^j)= \p{r^,z^,t-T)h{r^,z^,r^,z^,T)dr , 

—X 

where p is pressure, is the receiver range, is the receiver depth, is the source 
range, z^ is the source depth, and t is time. The Fourier transform of the equation gives 
the resulting pressure at the receiver in the frequency domain 

P{r,,z^,f) = P{r^,z^,f)H{r^,z^,r^,z^,f) < 2 ) 

where P and H are Fourier transforms of the functions in (1) such that 

X 

p{r^,z„t) = j P{r^,z^,f)e~^^”^' df p) 

—X 

and 



* The assumption of time invariance would at first seem to be unrealistic for the 
ocean. However, the ocean need only be quasi-time-invariant, i.e. the properties of the 
medium do not change over the time it takes to propagate from source to receiver. This is 
often referred to as the "frozen approximation". 
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The formulation of the ocean transfer function becomes critical and the best approach to 
this end remains in the solution of the Helmholtz equation. 

The Helmholtz equation governs the sound pressure field excited by a point 
harmonic source, therefore the Helmholtz equation evaluated at different frequencies will 
determine the transfer function. Assuming the complex pressure p(t,r) satisfies the 
pressure release boundary condition p(t,r)=0 at the surface and the outgoing radiation 
condition at infinity, the reduced wave equation in the water column becomes 



where p(t,r) is the acoustic pressure (Pascals) and K is the complex wave number. 
Assuming that the acoustic pressure has a time harmonic dependence, then substituting 
this acoustic pressure into the wave equation results in the time independent Helmholtz 
equation. The solution to this equation can be obtained by the parabolic equation method 
to approximate the solution. 

Historically, the parabolic equation method has been limited to modeling narrow 
angle sound energy propagation. More recently, parabolic equation acoustic algorithms 
have modeled sound propagation with propagation angles of up to 90 degrees from the 
horizontal. “>'2 Additionally, the technique readily implements models for coarse, rugged 
bathymetry and/or variable density stratification. In this thesis, the Collins' parabolic 
equation approximation method'^ computes the ocean transfer function. 
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The Helmholtz equation can be solved in cylindrical coordinates by removing the 
cylindrical separation effect terms and assuming negligible azimuthal variations. By 
defining O = -JT‘P and substituting into Eq. 5, the Helmholtz equation can be factored in 
out-going and incoming solutions. The out-going equation is 



^ = jk,^\ + XQ 
or 



( 6 ) 



where 







V 




1 dp d ^ 
p dz dz j 



( 7 ) 



This algebraically-difficult square root operator can be approximated by a series 
expansion. 

The Finite Element Parabolic Equation (FEPE)’ • of Collins uses a family of Pade 
series'"* to create higher order parabolic equations which are accurate for propagation 
angles close to the vertical. With Pade series expansion the square root operator becomes 



where 



vr+x-1 



h 1 + hn^ 



( 8 ) 






2 . T ifr 

sm“ 

2n + \ 2n + \ 



b^. 



= cos" 



IK 

2n + 1 



where n is the number of Pade terms. The number of Pade terms used determines the 
angle of sound propagation which can be accurately modeled, with excellent results up to 
90 degrees for //=4. 
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To derive the parabolic form of the complex pressure, O is expressed as 



Q=4i P= 



( 9 ) 



By solving for P and equating this to the ocean transfer function, Eq. 9 becomes 




( 10 ) 



Now substituting into Eq. 2, the frequency domain pressure signal becomes 




Taking the Fourier transform, combining exponential terms and defining t' = t — rj C q , 
the resulting received pressure signal becomes 



which can be solved efficiently with the inverse fast Fourier transform (IFFT). 

To solve the Helmholtz equation using the Collins FEPE, several assumptions are 
made. Specifically, the acoustic pressure has time-harmonic dependence as mentioned 
above. In addition, the sound source is an omnidirectional point source and backscatter is 
neglected. 
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B. THE OCEAN GRID 

The localization algorithm requires two assumptions with regard to the a prion 
knowledge of the source One, the source must have been propagated from a location in 
the ocean waveguide along a specifically known azimuth. The second assumption requires 
that the source emanate from a location on the superimposed grid. Figure 1 depicts 
representative grid. Considering that most common receiver is an array of hydrophones, 
the first assumption takes advantage of beamforming to localize in azimuth. Later 
discussion will expand on the second assumption requirement. These two assumptions 
lead to the creation of a 'grid' of possible source locations with axes of depth (z) and range 
(r) along the azimuth given. 

C. RECIPROCITY 

The final aspect to consider is that of reciprocity. With each grid point there is an 
associated transfer function from that point to the receiver. If it can be claimed that a 
signal sent from point a to point b in the ocean will appear at b the same as a signal sent 
from point b to a would appear at a, then reciprocity is consider to exist'^ Symbolically, 
this is 

or 

H{r,,z^,r^,z^,f)- H{r^,z^,r^,z^,f) ^ (, 4 ) 

where the H indicates the transfer function in the forward direction and the H indicates 
the transfer function in the reverse direction. This assumption is fairly accurate since the 
ocean and surrounding environment is for the most part reversible, except in the presence 
of strong currents. 
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Surface 




Figure 1: Example of grid and source-receiver relationship 
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III. LOCALIZATION METHODOLOGY 



The localizer works by using reciprocity and transmitting the original source signal 
back from the receiver to all the points in the grid. The resulting 'reverse' transmissions 
are then correlated with the original received signal, i.e., H{r,z)S is correlated with 

D= HS for all ("/* ,r) in the grid, such that 









( 15 ) 



The location with the maximum cross correlation will be the original source location. This 
method is dependent on knowing the original source signal characteristics. 

A variation in the algorithm results when lacking characteristics of the original 
signal. In this case, broadband signals can be localized by assuming that the transient 
signal will only result in a scaling of the transfer function in the frequency domain. The 
signal being an impulse in the time domain will be a constant in the frequency domain. 
Therefore, the transfer function H at all points in the grid are correlated with the 
received signal D = HS , 



Romir,Z\T) 






( 16 ) 



This will not result in a perfect correlation, but the maximum value of the correlation 
function at each point in the grid should reveal the source. Again, the location of the 
maximum cross correlation is the source location. If several relatively similar maximums 
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exist then multiple hydrophone data may be employed to discern the true source location^. 
This variation of the method will work for very short duration transients, on the order of 
100 msec or less. 
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IV. LOCALIZER APPLICATION 



A. MONTEREY BAY EXPERIMENT 

As the Collins FEPE is a prototype code, verifying its accuracy required a test to 
compare with a known signal transmitted in a known ocean waveguide. While solving the 
forward modeling problem, Westreich’^ tested the code in a shallow water acoustic 
tomography experiment in Monterey Bay. In December, 1988, a specially designed 
transmitter produced phase-encoded acoustic signals continuously for four days in the 
waters off Monterey, California to investigate the feasibility of using tomography for 
ocean interior and surface monitoring in a coastal environment.^"^ This environment 
presented problems since other methods of predicting the resulting signal at the end of the 
waveguide would breakdown due to the complexity. The FEPE code accurately modeled 
the time domain structure of the signal at the receiver when compared to the measured 
arrivals (Figures 4 and 5). Therefore, in comparison with the measured signal data, the 
FEPE modeled signal data agreement gave a high degree of confidence in the FEPE 
algorithm. 

The first series of computer experiments for this thesis investigated target 
localization using the Clay algorithm with the ocean environment of the Monterey Bay. 
The Collins FEPE code used the same range-independent sound speed profile and range 
dependent bathymetry as the Monterey Bay tomography experiment (Figures 6 and 7). 
The modified FEPE code modeled the reciprocal transmission of a Blackman pulse source 
from the Station J receiver location to predict time domain pressure signals at multiple 
depths and ranges. The localizer correlated these time domain pressure signals with the 
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Figure 4: Modeled received signal for Monterey Bay experiment.'^ 
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Figure 5: Actual measured signal at Station J in the Monterey Bay experiment. 
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previously modeled received signal for ranges of 30 to 55.75 km from the receiver on a 
line of bearing to the true source location. The correlator conducted checks on a grid of 
1 000 meter range increments by 1 5 meter depth increments. If the grid was appropriately 
chosen and assuming acoustic reciprocity holds, the algorithm should reveal the true 
source location at the position of the maximum value of the correlation. Figure 8 and 
Figure 9 show the maximum values of the normalized correlation coefficient plotted vs. 
range and depth. The correlation coefficient has its maximum value of 0.92 at the true 
source depth and range of 55.75 km and 885 m depth. The next highest correlation is 
0.76. The maximum value of this model-model correlation was not 1.0 because FEPE is 
only approximately reciprocal. The degree of reciprocity of FEPE and its sensitivity to 
different environmental scenarios is an important topic for future research The complex 
environment of Monterey Bay reveals a unique target location for a known source. 

The localizer results shown in Figure 8 and 9 also show a pronounced maximum. 
This leads to a problem with the localizer since it is very resolution dependent. As seen in 
Figure 9, locations within one grid spacing in depth have 0.92 to 0.74 drop in correlation. 
Therefore the localizer is very dependent on the grid and how it is overlaid upon the 
environment. Rovero*^ has explored methods that do not suffer as significantly from this 
criteria. 
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Range 

Figure 6: Monterey Bay bathymetry used in FEPE model. 



Monterey Bay Souncd Velocity Profile 




Sound Speed 

Figure 7: Monterey Bay Sound Velocity Profile used in FEPE model. 
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Figure 8: Source localization in range. 



Moximum Cross Correlation ocross Depths 




Figure 9: Source localization in depth. 
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B. 



BARENTS SEA EXPERIMENT 



Here, it is shown that the Clay localization algorithm s relatively insensitive to a 
priori knowledge of the actual signal. An experiment was conducted using the bathymetry 
and ocean conditions of a tomography experiment scheduled for the summer of 1992 in 
the Barents Sea. Figure 10 shows the range-independent sound speed profile used as input 
to FEPE and Figure 1 1 indicates the bathymetry used as input. The bat hy metrically simple 
region has depth changes from 163 meters to 320 meters, a 0.2 degree slope, over a 56 
km range. The algorithm using a modeled signal found the unique source location with a 
correlation value of 0.88. The next largest coefficient value was 0.76 (Figures. 12 and 13). 

As in the Monterey Bay experiment, Figures 12 and 13 reveal that the correct 
location is found, in part, because the grid was placed properly and one of the locations 
checked was on the grid. This may restrict the algorithm's utility unless a more robust 
method that resists this dependence is developed. Assuming this condition is met, another 
criteria that required investigation is that of source signal knowledge. 

To determine signal sensitivity the same environmental model was used but a 

/X 

Blackman pulse centered at 500 Hz +/-50 Hz was used for the source estimate S in (15) 
as shown in 






R (r ‘ r) = 






(17) 



instead of the actual signal. The only commonalty between the "true" and "modeled" 
signal is the bandwidth of the signals (Figure 14). 

This approximation of the source lowered the maximum value of the correlation 
coefficient to 0.54 at the correct location. However, the next largest value of the 
correlation coefficient was 0.43. As shown in Figures 15 and 16, this initial modification 
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Depth (m) 



of the Clay algorithm clearly defined the source location without ambiguity This shows 
that very limited information may suffice when attempting to localize a source in shallow 
water environments. 



Barents Sea Battny m e t ry 




Range (mxlO^) 

Figure 10: Bathymetry used for Barents Sea FEPE model. 
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Barents Sea Sound Velocity Profile 




Sound Velocity (rn/s) 



Figure 11: Sound velocity profile used for Barents Sea FEPE model. 




Range ( ) 

Figure 12: Localization in range for a known signal used when correlating. Source at 
1 5,000 meters range shown. 
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DeptH ( rnn) 

Figure 13: Localization in depth with a known signal used for correlation. 




Rrepueocy (1— Iz) 

Figure 14: Comparison between "real" signal spectral components and Blackman 
window model of signal used in Barents Sea computer experiment. 
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Figure 15: Localization in range with a source assumed as Blackman pulse in the 
frequency domain between 450 and 550 Hz. Source at 15,000 meters range shown. 




Depth! 

Figure 16: Localization in depth for a source assumed to be Blackman pulse in frequency 
domain between 450 and 550 Hz. Source depth 100 meters depth shown. 
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C. GENERIC EXPERIMENT 

With the knowledge that for short duration time domain signals, the source signal 
could be modeled as an impulse in time, a test case was generated to determine a lower 
bound on the necessity of source signal knowledge. The frequency band occupied by a 
signal transmitted through the linear time-invariant ocean waveguide will not change 
Therefore, a signal that has band limited characteristics could be correlated only over that 
band. In addition, the signal could be modeled as a constant amplitude over those 
frequencies. In these experiments, both a Clay localizer and a unknown source localizer 
are demonstrated for comparison 

1. Flat Waveguide 

The source is modeled as a Blackman pulse 224 Hz +/- 16 Hz in a 375 
meters deep, 10,000 meters long flat waveguide and with the sound velocity profile shown 
in Figure 17. The source was located at 150 meters and the receiver at 100 meters. The 
Clay localizer successfully located the source as expected, with a peak of l.O for the 
correlator value at 150 meters and 10,000 meters (Figures 18 and 19). The reciprocity of 
FEPE as demonstrated by the near perfect correlation is due to the range independence of 
the waveguide. The correlation with the impulse response of the ocean waveguide 
resulted in another accurate localization with suppressed peak of 0.692 (Figures 20 and 
21), which was still well above the side-lobe peaks. 
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Figure 17: Sound velocity profile used for generic waveguide FEPE model. 
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Figure 18: Localization in range for a known signal. Source at a range of 10,000 meters. 
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M o s< i pnn Lj nn of Cross — correlation Across Dopthis 




Figure 19: Localization in depth for a known signal. Source at a depth of 150 meters, 



Moximum of Cross — correlation Across Depths 




Depth 

Figure 20: Localization versus depth for an unknown signal (Source signal assumed as 
impulse in time). Source depth again shown to be 1 50 meters. 
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Figure 21: Localization in range for an unknown signal (Source signal assumed to be 
impulse in time.) Source shown to be from 10,000 meters range. 

2. Sloped Waveguide 

Figures 22 and 23 show the results of the localizer with no a priori source 
information but in a modified environment from above. In this case FEPE used the same 
sound velocity profile and a sloping bottom bathymetry, one meter rise for every 500 
meters in range, as input. Here the algorithm clearly defines the source location at 10,000 
meters range and 150 meters depth. The sharpness of the peak relative to surrounding 
depths and ranges makes apparent the requirement that the grid must map onto the source 
location for this algorithm to be successful. 
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Sloped bottom cose — maximum volue at each depth 




Figure 22: Sloped bottom case, maximum correlation for each depth strata with source 
depth resolution shown. 



Sloped bottom case — maximum value at each range 




Range 

Figure 23: Sloped bottom case, maximum correlation values for each range increment. 
Source location resolved at 10,000 meters range. 



25 



3. Rugged Waveguide 

The final case examined used a more rugged bottom profile as shown in 
Figure 24, This case also resulted in accurate source localization without knowledge of the 
specific source spectrum. Figures 25 and 26 show the source location clearly defined at 
1 50 meters depth and 1 0,000 meters range. These do show, however, unexpectedly high 
values for the localizer for points other than the original source location. Since this 
environment is more complex than in the prior cases, expectation was that the relative 
transfer functions would resolve the actual location clearly. 
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Figure 24: Rugged bathymetry case used in FEPE model input for Generic waveguide 
tests. 
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Unknown source locolizer 




Figure 25: Source localization for rugged bathymetry with no source knowledge 
maximum values at each depth strata. Source revealed at 150 meters depth. 



Unknown source localizer 




Ron g e C 

Figure 26: Source localization for rugged bathymetry with no source knowledge 
maximum values at each range step. Source appears to be at 10,000 meters. 
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V. 



CONCLUSIONS 



This thesis has shown by computer simulation that transient signals in shallow 
water environments can be eflfectively localized. Scenarios used included Monterey Bay, 
Barents Sea, and generic waveguide simulations in which the signal source was both 
known and unknown 

Although these were not entirely new methods, this thesis has demonstrated that 
the short duration time-domain transient localizers work with the Finite Element Parabolic 
Equation model of Collins. With a well known environment, the algorithms shown have 
clearly and consistently revealed the transient source location for a variety of cases 
without ambiguity. Moreover, these investigations demonstrated that source signal 
parameters need not be completely specified for successful localization. Critical to these 
developments were the especially important new results that came from the use of the 
FEPE model to develop the ocean transfer function in the shallow water environments 
studied. 

Future research should include determination of optimum FEPE grid requirements 
to accurately model transmitted signals when compared to experimentally measured 
signals. Also to achieve tactical Navy utility*^, the ocean modeling routines and 
localization algorithms require additional investigations into the effects of environmental 
uncertainty^*’ and multiple hydrophone data sets. Computational demands should be 
explored to determine the best computer architectures for executing these algorithms. 
Finally, additional topics for future research should include investigation of source 
localization in noisy environments, integration with detection and classification algorithms, 
and an investigation of the performance bounds of the FEPE reciprocity. 
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